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Abstract: Using the full radiation transfer function, we numerically calculate the CMB 
angular bispectrum seeded by the compensated magnetic scalar density mode. We find 
that, for the string inspired primordial magnetic fields characterized by index ub = —2.9 
and mean-field amplitude B\ = 9 nG, the angular bispectrum is dominated by two primor- 
dial magnetic shapes. The first magnetic shape looks similar to the one from local-type 
primordial curvature perturbations, so both the amplitude and profile of the Komatsu- 
Spergel estimator (reduced bispectrum) seeded by this shape are almost the same as those 
of the primary CMB anisotropics. However, for different parameter sets (Zi, this "local- 
type" reduced bispectrum oscillates around different asymptotic values in the high-^a regime 
because of the effect of the Lorentz force, which is exerted by the primordial magnetic fields 
on the charged baryons. This feature is different from the standard case where all modes 
approach to zero asymptotically in the high-Z limit. On the other hand, the second mag- 
netic shape appears only in the primordial magnetic field model. The amplitude of the 
Komatsu-Spergel estimator sourced by the second shape diverges in the low-/ regime be- 
cause of the negative slope of shape. In the high-/ regime, this amplitude is approximately 
equal to that of the first estimator, but with a reversal phase. 
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1. Introduction 



In the inflationary scenario Q, the quantum fluctuations of the scalar field(s) are re- 
sponsible to generate the initial conditions for the Cosmic Microwave Background (CMB) 
anisotropies. The current observations Q from large scale structures are consistent with 
an almost scale invariant, Gaussian primordial density perturbations generated during in- 
flation. However, with the improvements of measurement precision, any small deviations 
from the Gaussian distribution enable us to distinguish different cosmological models. Like 
the role colliders play in particle physics, measurements of non-Gaussian features provide 
microscopic information on the interactions of the inflatons and/or curvatons. Constrain- 
ing and detecting non-Gaussianity (NG) have become one of the major efforts in modern 
cosmology. A variety of potentially detectable forms of primordial non-Gaussian features 
from inflation models have been intensively investigated (see |^, |^ for a review). The ef- 
fects of primordial non-Gaussian curvature perturbations on CMB anisotropies have also 
been studied in the recent papers H, 0, |^ 10, 11, 12, [l^, 0, |l5[. The current lim- 
itations on the primordial bispectra from WMAP-7yr data are — 10 < /nl**' < 74 and 
—214 < /n'l''' < 266 at 95% CL ||2|, where /^l*^^ and /^l"'' are the non-hnear parameters 
of the "squeezed" and "equilateral" momentum configurations, respectively. 

Except for the possible NG from the inflationary dynamics, the primordial NG might 
come from other mechanisms. One interesting possibility is that the non-Gaussianities are 
sourced by the primordial magnetic fields (PMFs) in the large scale structures ||l^, 17, 18|. 
The astrophysical observations about the spiral/elliptical galaxies and rich clusters indicate 
that our universe is permeated with large scale coherent magnetic fields with the magni- 
tudes ranging from hundreds of nG to few fiG |1£, 2C, 21, however, their origins are 
still not yet fully understood. The dynamo mechanism explains the origin of the galactic 
magnetic fields with amplification of a small frozen-in seed field to the observed /xG field 
through turbulence and differential rotation 1 23 , ^] . And the gravitational adiabatic com- 
pression may generate the magnetic fields in clusters during the collapse of a protogalactic 
cloud [25, 26, 27]. Cosmological phase transitions in the early universe may produce the 
tiny magnetic seed fields, which are required by the above mentioned amplification mecha- 
nisms, such as the electroweak phase transition |28, QCD phase transition |^, ^ and 
the inflation with the broken conformal invariance [32|. 

In recent years, intensive effort has been devoted to studying the imprints of magnetic 
fields on the CMB anisotropies, which are nicely reviewed in |3^. The contributions to 
the CMB angular power spectrum from the scalar perturbations induced by PMFs are 
investigated in |||, ||, H, |3|, |o|, |l], |4|, ||, ||, ||, |4^, |8| , from the vector perturbations 
in 50, 51, 43, 44| and from tensor mode in [51, 52, 53, 4^, 44], respectively. Some other 
phenomena induced by PMFs, such as Faraday rotation, damping of Aflven waves, effects 
of PMFs on seeds for large scale structures and on neutrino masses are investigated in 
]54], [55, 56], ]57] and ]58, 59], respectively. And the new constraints on PMFs from CMB 



anisotropy and large scale structure data are reported in ]6C, 61 



In the inflationary scenario, the NG signals come from the high order curvature per- 
turbations. However, even at the lowest order, PMFs can still generate some non-Gaussian 
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features in the CMB anisotropies, since the magnetic energy density and anisotropic stress 
induced by PMFs are naturally non-Gaussian variables. Such signatures have been investi- 
gated in 1 63, 64, 35, 66, |5^, 69], but for the homogeneous magnetic fields with fixed 



direction which break the spatial isotropy and result in the north-south asymmetry on the 



CMB sky. However, as pointed out in |17, the stochastic PMFs are able to generate 
a distinctive non-Gaussian signal in the CMB anisotropies with an amplitude comparable 
with the one from the primary curvature perturbations. 



The authors in |17, analytically calculate the CMB bispectrum from the stochastic 
PMFs, but only in the Sachs- Wolfe regime {I < 10). In this paper we calculate the angular 
bispectrum from scalar perturbations induced by PMFs with the full transfer function. 
Because the scale-invariant magnetic power spectra are strongly inspired by string cosmo- 
logical model as a consequence of the breaking of conformal invariance during the pre-big 
bang phase |3^, 37], following most of the literatures in this subject (for example, — — 3 
in [|^] and ub = —2.9, ±2 in [|l8|), in this paper we take the magnetic index ub = —2.9 
and mean-field amplitude B\ = 9 nC. In this model, we find that the angular bispectrum 
is dominated by two primordial magnetic shapes. The first magnetic shape f^^\k,q,p) 
( |5.9| ) looks similar to the one from local-type primordial curvature perturbations, so both 
the amplitude and profile of the Komatsu-Spergel estimator (reduced bispectrum) seeded 
by this shape are almost the same as those of the primary CMB anisotropies ] ]70| , |105| , 
(see Figure 0, ^, and ^). However, for different parameter sets {h,h), this "local-type" 
estimator b^^j^^.^ ( |5.18| ) oscillates around different asymptotic values in the high-Zs regime 
because of the effect of the Lorentz force, (see Figure ^ and This feature is different 
from the standard case where all modes approach to zero asymptotically in the high-Z limit. 
On the other hand, the second magnetic shape f^'^\k,q,p) ( ^.10 ) appears only in the pri- 



mordial magnetic field model. However, the amplitude of the Komatsu-Spergel estimator 
^[li^s (EjH) sourced by the shape f^'^\k,q,p) diverges in the low-/ regime because of the 
negative slope of shape. In the high-/ regime, this amplitude is approximately equal to 
that of the first estimator b^^j^h^ with a reversal phase, (see Figure ^, 10 and 11). 



The rest of this paper is organized as follows. In section |2|, we firstly present the 
Maxwell and conservation equations which govern the behaviors of electromagnetic fields in 
the curved spacetime. Then we calculate the primordial magnetic power spectrum induced 
by PMFs under the ideal magnetohydro dynamics approximation. The linearized scalar 
equations for each individual matter component in the Cold Dark Matter (CDM) model 
and the gravitational fields are given in section ^. In section ^, we derive two magnetic 
initial conditions in the deep radiation dominant era, and then calculate the CMB angular 
power spectrum numerically by using these initial conditions. The numerical calculations 
about CMB bispectrum signatures seeded by the compensated magnetic density mode are 
analyzed in section ^. Finally, we conclude in section |6[ 

2. Stochastic primordial magnetic fields 

In this section, we firstly present the Maxwell and conservation equations which govern the 
evolution of electromagnetic fields in a curved spacetime. Then we calculate the primor- 
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dial magnetic power spectrum induced by PMFs under the ideal magnetohydrodynamics 
approximation . 



2.1 Electromagnetic field in a curved spacetime 

In this subsection, we present the Maxwell's equations and conservation equations in a 
covariant formulism^. Following the formulism, the electromagnetic (Faraday) tensor Fab 
can be decomposed into an electric and a magnetic component as 

Fab = 2u[aEb] + SabcB", (2.1) 

where Ea = FabU^ and Ba = EabcF^^ /2 are respectively the electric and magnetic fields 
experienced by the observer with 4- velocity u"" {Eau'^ = BauP" = 0). The Faraday tensor 
also determines the energy-momentum tensor of the electromagnetic field as^ 



T. 



(em) 



ab 



1 

47r 



-FacF% — -^FcdF^gab 



(2.2) 



Combining ( |2.lD with ( |2.2D , we arrive at the irreducible form of t'^'^ 

1(^2 ^ B^)UaUb + \{E^ + B^)hab + 2q^aUb) 



rp{em) _ ^_ 

«6 ~ 4vr 



^If • (2-3) 



Here E'^ = EaE°- and B^ = BaB"^ are the square magnitudes, qa = EabcE^B'^ and vr^ 
{—E(^aEb) — B(^aBb))/^'^ are the electromagnetic Poynting vector and anisotropic stress 
tensor, respectively ^ . In this paper, the round, squared and angled brackets denote the 
symmetric, anti-symmetric, and symmetric trace-free parts of a tensor, respectively. 
In the standard tensor form the Maxwell equations read 

VbE^"" = r , V[,F,,] = ^ r?"'^''Fb,;rf = , (2.4) 

where J" is the 4-current that sources the electromagnetic field. With respect to the 
Ua-congruence, the 4-current splits into its irreducible parts according to 

r = nu"" + J" , (2.5) 

with n = — Jau"', J'^ = h°-f^-fi and Jau"" = 0. By virtue of the irreducible form of Fab and 
J"", the timelike parts of the Maxwell equations read 

ha'E, = -\QEa + {Oab + Eabc^lE'' + EabcA^B' + CUrl Ba - Ja , (2.6) 

ha'B, = -\QBa + {aab + eabc^'')B^ " ^abc^^^' " curl Ea , (2.7) 



^The covariant approach to cosmological perturbations is shortly reviewed in Appendix ^ and the 
definitions of covariant variables, such as 4- velocity u", expansion rate O, etc. can be found there. 
^In this paper we take the unit conventions as c = ft = Mpi = l/SvrG = 1. 

^Comparing with the convensional definition about the electromagnetic anisotropic tensor tt^^', such as 
the one in the Jackson's textbook ||7^ , the definition in this paper is different from the convensional one by 
a factor of l/47r. 
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while their spacehke components provide the constraints 



DaE^ + lUaB" = fl , (2.8) 
DaB" - loOaE" = , (2.9) 

where Da denotes the spatial derivatives with respect to projected metric hab and its 
definition is presented in ( |A.2| ). 

Besides, the 4-current conservation law S/aJ"" = gives the continuity equation of 
charge density 

A = -G/X - DaJ"" - AaJ"" . (2.10) 

The equations (p.6D,( ^^ ,( |2.8| ),( |2.9| ),( |2.1[l| ) form a complete set of equations which evolve 
the electromagnetic field in a curved spacetime. 

2.2 Ideal MHD approximation in the Universe 

A good conductor throughout the history of the Universe allows us to study the electro- 
magnetic field in the universe within the limits of ideal magnetohydrodynamics (MHD) 
approximation. By means of Ohm's law, the spatial currents Ja read 

Ja = , (2.11) 

where ? represents the scalar conductivity of the medium. The MHD approximation states 
that, in the limit ? — )■ oo, we can neglect the electric field Ea- Hence, the energy-momentum 
tensor of the residual magnetic field becomes 



^B'^UaUb + ^B'^hab 



+ ^ib\ (2-12) 



with the anisotropic tensor vr^^^ = —Bi^^Bb) /4vr. From the above expression we can identify 
the energy density of PMFs as A^^^ = B'^/Stt. In addition, the Maxwell equations reduce 
into a single propagation equation 

Bta) = i^ab + eabcto" " lohab] B^ , (2.13) 



3 



and three constraints 



Ja = curl Ba + EabcA^B" , (2.14) 

/i = 2uj''Ba , (2.15) 
= D^Ba . (2.16) 

2.3 Primordial povi^er spectrum induced by PMFs 

In this subsection, we present the primordial power spectrum induced by scalar perturba- 
tions from the stochastic PMFs. In the local rest frame u"" = (1,0), Bau"' = leads to a 
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vanishing temporal component of Ba and then we have 



1 

47r 



B\t,iL)Bk{t,^)5ij - Biit,x)Bjit,x) 



(2.17) 

(2.18) 

(2.19) 



Furthermore, in the ideal MHD regime we can separate out the time evolution of PMFs, 
Bi{t, x) = Bi{-x)/a?. Hence, in what follows we concentrate on the time independent spatial 
component -Bi(x) and take them as statistically homogeneous and isotropic random fields. 
The transversal nature of PMFs leads to 



■ Pi 



Si(k)m(k')) = (27r)3^p(^)(A;)5(k - k') , k < kn , 



(2.20) 



where Py = 5ij — kikj is the projector onto the transverse plane, ko is the wavenumber of 
damping scale and p(^) (A;) is the primordial magnetic power spectrum. For some specific 
magnetogenesis models P^^^ (k) takes the power law form 



P^^\k) = Ak'''' . 

In the above expression, we have adopted the Fourier transform convention as 



A(x) 



1 



d^kBi{]s.)e 



— ikx 



(27r)3 

Bi{k) = J d3^5i(x)e*-^ 



(2.21) 

(2.22) 
(2.23) 



It is convenient to introduce the Fourier components of the PMF energy density contrast 
A^^^ and scalar part of the anisotropic stress tensor tt^^ as 



(B)„-ik-x 



TT 



{B)i 



(x) 



(27r) 



d^k vrP) 



-6' J - k'kj 1 e- 



jk-x 



(2.24) 
(2.25) 



Thus, we obtain the expressions for A^^^ and tt^^ from the momentum convolution 



TT 



8tt J (27r)3 
d^p 



B\p)Bi{k-p) , 



Sit J (27r)3 [ 



1 



hB\p)k^Bj{\^ - p) - -B\p)Bi{\^ - p) 



(2.26) 
(2.27) 



Since we are interested in the PMFs in the linear perturbation regime, we therefore define 
the magnetic comoving mean-field amplitude by smoothing over a Gaussian sphere of the 
comoving radius A = 1 Mpc (/^ = e"^'^^^/'^) as 



(i?,(x)i?,(x))U = i?f . 



(2.28) 
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For the power law model (|2^ ), can be given by the Fourier transform of the product 
of the power spectrum P^^\k) and the square of the filter transform fk, 



Bl 



(27r)3 



d'kP^^\k)\h\ 



2A 



1 ^ (uB + 'i 



(2.29) 



(27r)2 A"s+3 

where we require the spectral index ub > —3 to prevent the infrared divergence at the 
power spectrum level. Plugging ( |2.29D into (2.20), we arrive at 



3P,, (27r)"i^+5^f fena 



(S.(k)B*(k')) = (27r)=^ 



5(k - k') 



k < , 



(2.30) 



2 ; "-A 

where fc^ = 27r/A. For all scales smaller than the damping scale {k > kn — 4.5 Mpc^^) the 
spectrum vanishes. 

Furthermore, we can obtain the two-point correlation functions for Aj, and vr^ by 
using the Wick theorem 

(A(^)(k)A(^>(k')) = ^^^^ / d^p(^\p)p(^H\k-p\){l + ^^') , (2.31) 
(vr(^)(k)vr(^)*(k')) = ^^1^ / d'p P(^)(p)p(^)(|k-p|) 

^ " (2.32) 



^ o 1 2 

27P/^+ 4/^ 



where ^ = p-(k — p), j = k ■ p and /3 = A; • (k — p). Following most of the literatures in 
this subject, and also due to that the nearly scale-invariant magnetic power spectra are 
strongly inspired by the string cosmological models as a consequence of the breaking of 
conformal invariance during the pre-big bang phase 37 1, we take ub = —2.9 in the 



following calculations. Ignoring the cutoff in the definitions of p(^) allows us to integrate 
( |23ll) and (|23^ ) semi-analytically ||] 

2 



P-KB (^) 



42.37 
16 

9 X 14.55 



(27r)"fl+252 



2r(^^)p(7) 



(27r)"s+252 

2r(^^)p(7) 



^ \ 2ns +6 
2 ^ X 2ns+6 



where we have used the convention about dimensionless power spectrum 

{X{k)X*{k')) = 27r2(2^)35(k - k')k-''Px{k) , X = Ab, ttb 



(2.33) 
(2.34) 

(2.35) 



3. Basic equations 

In this section we present the magnetic linearized scalar equations in Fourier space^. Firstly, 
we define the scalar-valued harmonic function on the exact Friedmann-Robertson- Walker 
(FRW) background 

a^D^Q^^\k) + k'^Q^^\k) = , Q(0)(A;) = 0, (3.1) 

*The set of the linear equations for all matter components in the coordinate space can be found in 
Appendix @ 
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where a is the scale factor and the superscript (0) represents the scalar mode. The covariant 
temporal and spatial derivatives are defined in ( |A.2 ). Arming with the scalar harmonics, 
we can calculate the rank-/ Projected Symmetric and Trace-Free (PSTF) tensors by virtue 
of the recursion relation 

Q^>) = -lD^aMt.)^k). (3.2) 



Another useful relation is 

•,(0) 



(0)... k I 



QZS'') ' (3-3) 



a {21 - 1) 

where the constant K is related to the spatial geometry of the universe {K = 0, 1 
corresponds to a flat, closed and open universe, respectively). In the above expressions, we 
have used the covariant spherical multipole expansion 

oo 

f{x\p^) = Y,FAAx\E)e^^ =F{E) + Fa{E)e'' + Fai,{E)e^e'' + -- - , (3.4) 
z=o 

where the PSTF tensor reads Fai{E) = F(^aia2---ai)iE)- 

Next, we expand all dynamical variables in terms of the harmonic tensors, which is 
similar to the Fourier series expansion. For the multipoles of intensity brightness of photon 
Iai and neutrino Gai, we have 

.(0) T(o)ri.^n(o)^ 



k \n=0 J 

GA=GY.iU-^n^] sPmf^ik), 1>1. (3.6) 

k \n=0 / 

where / = p^'^\ G = p^^^ and k>^^ = [l — (/^ — 1 — 'm)K/k'^] ^^"^ for / > m. For the sake of 
briefness, we will suppress the scalar superscript (0) and momentum k implicitly in QA^^k), 
Il^\k) and Gi^\k) in the rest part of our paper. And for other gauge-invariant variables, 
we have 

= T = (3.7) 

k 

#=P^^^E'^?^2- (3-8) 

k 

y'i' = Y.^Qa, (3.9) 

k 

^» =p»^vr«Q„,, (3.10) 



k 

a 



Za = - Y^—ZkQa, (3.11) 

k 

-^£kQab ! (3.12) 
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E-O-fc Qab 
n 



^ — ' n 



ah 



4ri) 



(7) 



(3.13) 

(3.14) 
(3.15) 

(3.16) 



where the superscript (i) denotes photon (7), massless neutrino {y\ baryon (6) and Cold 
Dark Matter (c), respectively. Since p^^^ vanishes at the background level, in ( p. 16 ) and 
(3.15) we therefore normalize vr^^^ and A^^^ by photon density p^'~^\ In the above conven- 
tions we have 



(3.17) 
(3.18) 



We are now ready to derive the scalar multipole equations for all matter components. 



3.1 Photons 

From (|B.9| ), the complete Boltzmann hierarchies for the total intensity of photon are 

k 



+ 1) ^(0) -J I (0)^ 



{21 + 1) 



-UeCTT 



{21 + 1) 

10 



6 a 15 a 



'^i - IqSio - -^vl '611 - —12812 



(3.19) 



where h = {kZ^/a — QA^) /3, the dot is derivative with respect to cosmic time t and 
rieCTT is the differential optical depth of the Thompson scattering. The first three hierarchy 
equations are 



(3.20) 



for the monopole case (Z = 0), 



(3.21) 



for the diploe case {1 = 1), and 



. C-yI 3 A; ^ 2 k (y) 8 k 9 (-y) 

TT^^^ + --I3- --ql'' -—-ak = -— neCTTvr^^^ 
^ 5 a 5a ^ 15 a 10 " 



(3.22) 



for the quadrupole case, respectively. 
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3.2 Massless Neutrinos 



Because the massless neutrino only gravitate, the Boltzmann hierarchies for the total in- 
tensity is similar with the one for photons except that in the right hand side, the Thompson 
scattering term vanishes: 



Gi + 



+ 1) , (0) 



k 



a [{21 + 



+1 



k) 'Gi-i 



{21 + 1 



+4:hSio + ^-AkSn - ^-i^2^(Tk5n = 
6 a 15 a 



(3.23) 



Because the massless neutrinos behave like coUisionless relativistic particles, we treat them 
as the improved fluid, i.e. we need expand in the multipole series to octupole at least. So 
in what follows we list the first four hierarchies: 
monopole {I = 0) 



a V 3 



rOAfc = 



dipole (1 = 1) 
quadrupole (Z = 2) 

octupole (/ = 3) 
3.3 Bayrons 



4^^ + ^^(2-r-Ar+4^.) = 0, 



• M , 3fc 2k iy) 8 k 



GS = -^T^k 



k3 (u) 
a 7 



(3.24) 
(3.25) 

(3.26) 
(3.27) 



For baryons and CDM we use the fluid approximation and neglect their anisotropic stress 
tensors, i.e. we characterize baryons and CDM only by the energy densities and velocities. 
For baryon density contrast we have 



P 



k 



-(Zfe + i;f)-eA 



+ c; 



2 P 



(by 



eAf = 



(3.28) 



where we use Dap'^^^ = c^Dap^'^y The baryon velocity equation reads 



1 + 



P 



(by 



lib) 



rieCTTP 



a " 



(3.29) 



where the first term in the right hand side denotes the usual Thompson scattering and the 
second new term for the Lorentz force from PMFs. 
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3.4 Cold Dark Matter 

For CDM we have 

and 



(3.30) 
(3.31) 



respectively. 

3.5 Gravitational equations 



The evolution equations of gravitational field read 



Zu H Zi, H — — 

" 3 ^ 2k 



3a 
'2k 



Ak --Ak = 
a 



(3.32) 



3^'=+2fc L 



(p + p)cTfc + /.Wg« 



3 -p' 



(0 



TT 



2ifc2 







And the corresponding constraint equations are 



g(Z,-a,)+p«,f = 



(3.33) 
(3.34) 

(3.35) 
(3.36) 



4. Compensated magnetic initial conditions and CMB power spectrum 

In this section, we analytically extract the scalar modes seeded by PMFs in the deep 
radiation dominated era. Then we numerically calculate the CMB angular power spectrum 
by using two compensated magnetic initial conditions. 

4.1 Equations in the tight-coupling approximation 

In this subsection we propagate the covariant equations in the zero-acceleration frame, in 
which the CDM velocity vanishes. For the density contrast of different species, we have 



Ai^^' + ^(^^. + 4'^)=0, 
Ar + ^Q^. + 4^^)=0, 



(4.1) 
(4.2) 
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Af +/c(z, + .f)=0, 
A'j^^' + kZk = 0, 



(4.3) 
(4.4) 



where ' = d/drj is the derivative with respect to the conformal time rj. For simplicity, we 
also set the pressure and sound-speed of baryon fluid to zero {p^^^ = = 0) . 

In the deep radiation dominant era, photons are tightly coupled with ionized baryons 
through the Thomspon scattering. This allows us to deal with them as a single baryon- 
photon fluid with a common fluid velocity q^^^ . Furthermore, PMFs also exert the Lorentz 
force onto the baryon-photon fluid. Hence, under the tight-coupling approximation {q^^^ — 
q^) ~ Av^^ 1 3) , the velocity equation of the baryon-photon fluid takes the following form 



a(7) , 

3(1 + ^ 



Af ) = 0, 



(4.5) 



{\^Ry^ 3(1 + i?) ^ 3(1 + i?) 

where H = a' /a is the conformal Hubble constant and R = Ap^^^ /3p^'^^ is the photon to 
baryon ratio. 

Since the massless neutrinos behave as coUisionless relativistic particles, they can pre- 
serve the non-vanishing octupole signals 



A 



n^^y + lkgs-lkq^:^-l^ka,=0, 
The gravitational field equations in the zero-acceleration frame read 



2 r 



(B) 



.(b) 



(4.6) 

(4.7) 
(4.8) 

, (4.9) 



2k 



+ 



0, 



a'j. + UcTk + k£k + 



m 

~2k 



2 , 







(4.10) 



(4.11) 



where we define the density fraction as R^ = P'y/ Ru = Pv/ P-, RbV = Pb/p and RcT] = Pc/ P- 
Note that in our definitions i?^ and Rc have the dimension (length)"^. In addition, the 
gravitational constraint equations are 



28, 



9n^ 



i?,7rf ) + R^ttI^^ + R,A^^^ + R,Af^ + R^A^^^ + i?eAf + ii^Af 



R,q^^^ + R.iU RbvP 



= 0, 



-2fe — CTfe -I- 



R^q^^^ + RJ^^ + RbV^^^ 



= . 



(4.12) 
(4.13) 
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4.2 Compensated magnetic initial conditions 

Since the radiation species dominate our universe during the initial era {R-y + Ry ~ 1), 
usually one neglects the matter contributions when derives the adiabatic initial conditions. 
However, as demonstrated in |45, 46], one cannot neglect the matter contributions in the 
case of existence of PMFs due to the compensation mechanism between the radiation den- 
sity perturbations and those of PMFs. So it is essential to take the matter contributions 
into account when we derive the magnetic initial conditions. In addition, it turns out con- 
venient to introduce a new characteristic length scale Rm = Rb + Rc ~ Pmivo)/ \/ -ipr ivo) — 
5 X 10~^Mpc~^. In what follows, we list two different compensated magnetic modes in- 
cluding matter contributions. 

The density sourced compensated magnetic mode: 
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And the anisotropic stress tensor ir), sourced mode: 
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2„2 



krj , 



7 2 2 
-k 7] 



(4.29) 
(4.30) 
(4.31) 



84(15 + AR^ 

where rjs = —{28k + ^''fc/^) is the curvature perturbation in the synchronous gauge. 
4.3 CMB power spectrum 

By virtue of the above initial conditions, we formally integrate the set of evolution equations 



over the line of sight 0, @, g9|, |Tog, |ToT| 



4 / dte-^ 



k 3 , (o)s_i. 



-O-fc + T^rieCTT(4 ) ^2 
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-Ak - UeCTTVk -7-3l{x) 

a I ax 



1 fk 

-i-Zk- OAk 
6 \ a 



1 

-UeaTM 



(4.32) 



where t = J neCTdt is the optical depth, x = kx with x is the comoving radial distance 
along the line of sight and ji{x) are the spherical Bessel functions. 
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Figure 1: The CMB spectrum of TT mode with the magnetic index ns = —2.9 and comoving 
magnetic mean- field amplitude Bx = 9 nG. The red solid curve stands for the primary adiabatic 
mode, the green dashed one for the magnetic density sourced mode and the blue dotted one 

for the magnetic anisotropic stress tt^ sourced mode, respectively. 



Then we expand the temperature contrast {6t = ST/Tq) in the multipole series 

oo oo 

^Tie") = j^Al'lA,e^' =7rY,Y.^i'^^STi{k)QA,e^' , (4.33) 
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Figure 2: The CMB spectrum of EE mode with magnetic index ns — —2.9 and comoving magnetic 
mean-field amphtude B\ = 9 nG. The red soHd curve stands for the primary adiabatic mode, the 
green dashed one for the magnetic density sourced mode and the blue dotted one for the 

magnetic anisotropic stress tt^, sourced mode, respectively. 

where in the second equality we rewrite the multipole coefficient Xi{k) = CkgTi{k) with 
the transfer function gTi{k) and random variables C^, which source the CMB anisotropies 
with primordial power spectrum 



{CkCl,) = C\k)5kk' 



(4.34) 



Usually Cfc are the primordial curvature perturbations, however, in this paper they are the 
density contrast or anisotropic stress vr^ of PMFs, and their power spectra are given 
in p:33D and (|^). 

Armed with the primordial power spectrum, we finally obtain the CMB angular power 
spectrum 



Ci 



IT 



dink C\k)\gTiik)f 



(4.35) 



In Figure and H, we plot the CMB TT, EE and TE spectra with the primary adi- 

(B) 

abatic mode (red solid curve), magnetic density sourced mode (green dashed one) 

( Fi) 

and magnetic anisotropic stress vr^ sourced mode (blue dotted one), respectively. In our 
numerical calculations, we modify CAMB code |102[| and set the amplitudes of primordial 
curvature perturbations to unit. The results show that PMFs contribute a tiny part to the 
CMB power spectra, however, in the next section, we will demonstrate that they will give 
a dominant contribution at the bispectrum level. 

5. CMB bispectrum 

In this section we numerically calculate the CMB bispectrum seeded by the compensated 
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Figure 3: The CMB spectrum of TE mode with magnetic index ub = —2.9 and comoving magnetic 
mean-field ampUtude B\ = 9 nG. The red solid curve stands for the primary adiabatic mode, the 
green dashed one for the magnetic density sourced mode and the blue dotted one for the 

( B) 

magnetic anisotropic stress tTj, sourced mode, respectively. 



magnetic density mode. 
5.1 Analytic formulas 

Firstly, let us shortly review the analytic formulas to calculate CMB angular bispectrum 
1 105 , 104 , [fo| . In ( [4.33 ) we decompose the temperature contrast St by the covariant 
approach, now we expand it in terms of the spherical harmonics, which are more familiar 
to us, 

'^r(n) = ^ ai„iYim{n) , (5.1) 

Im 

where n denotes the unit direction vector. The CMB angular bispectrum is defined as 

where B^f^^^^ must satisfy the triangle conditions and the selection rules: mi+m2+m^ = 
0, ^1 + ^2 + ^3= even and \li — lj\ < Ik ^ h + Ij for all permutations of indices. Note that 
Gaunt integral satisfies all the conditions mentioned above 



(2Zi + l)(2/2 + l)(2/3 + l) Ml /2 ^3 \ hhl 



3 



(5.3) 



47r \ / \ nil rn2 I 

where the matrices denote the Wigner-3j symbol. Therefore it is convenient to introduce the 
reduced bispectrum (Komatsu-Spergel estimator) bi^i^i.^ 11^1 to replace B^f^^"^^ without 
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any loss of information 
Thus, the observable angle-averaged bispectrum can be written as 



KLTr' = Guirr'Ki.i^ ■ (5.4) 



S = y I h h h \ ^rn,n^,m, 

'1*2(3 \mim2m3 'l'2«3 

mim2m3 \ ^ ^ / 

(2/i + l)(2/2 + l)(2/3 + l) / h h h 
47r I 



(5.5) 



In order to calculate the reduced bispectrum h-^i^i^, we need obtain the form of primor- 
dial bispectrum F(ki,k2,k3)- For a slow roll inflation model, the non-gaussian curvature 
perturbations are usually parameterized by a single constant parameter and the 
Gaussian random variable Cl(x) the real space 

C(x) = a(x) + /]vT'[cl(x) - (C'(x))] . (5.6) 

In the Fourier space, the local-type primordial curvature bispectrum (("(ki)("(k2)C(k3)) = 
F(^{ki, k2, fe3)(5(ki + k2 + ks) can be obtained by performing the momentum convolution 

F^ikiMM « fNt{Pd^i)Pd^2)kl + {ki,k2, ks) perm.} , (5.7) 

where is the primordial power spectrum of curvature perturbations (C)- However, we are 
interested in the PMF signals in CMB bispectrum which have the essential non-Gaussian 
characters. The magnetic density contrast bispectrum (A^^^(ki)A('^^(k2)A(^^(k3)) = 
{ki, k2, A;3)(5(ki + k2 + ks) has been derived analytically in [|^, |l^ 

FAAk,q,p) = ^g^2p(7)3 [f^'Hk,q,p) + f^^\k,q,p) + f^^\k,q,p)\ , (5.8) 
f^'Hk,q,p) = ""f fc^"^+V"+V + {k,q,p) perm. , (5.9) 

f^^\k,q,p) = „^ , k^q^p^ + {k,q,p) perm. . (5.11) 
Sub + o 

From the above expressions, we notice that the first magnetic shape f^^\k,q,p) looks 
similar to the one from local- type primordial curvature perturbations ( |5.7D , if the magnetic 
index takes the nearly scale-invariant value ub ~ —3. So, we will abuse the phrase "local- 
type" for the magnetic shape f^^\k,q,p) in this paper. 

Now, we are ready to calculate the CMB bispectrum sourced by the magnetic density 
contrast A(^)(k). Following the standard procedure |103|, the reduced bispectrum can be 
expressed as 



hihh = (8^)^ / 2;^(ix / dink / ding / dlnp ji^{kx)ji^{qx)ji^{px) 
Jo Jo Jo Jo 

y.F^g{k,q,p)gTi^{k)gTi^{q)gTi^{p) , (5.12) 
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where gTi{k) is the transfer function and x is the comoving radial distance along the line 
of sight. From (pD-(pID, we can see that the integral ( |5.12[ ) is determined by four kinds 
of momentum integrations 

^ t'' dink k'^^^^3i{kx)gTi{k) , (5.13) 
Jo 

hf\x) = f " dink k''^+^ji{kx)gTi{k) , (5.14) 
Jo 

^(7)r^^ _ / 



bY\x)= dink k'jiikx)gTi{k) , (5.15) 







bf\x) = [ dink k^''^^^ ji{kx)gTi{k) . (5.16) 



Then, we can express the reduced bispectrum ( |5.12 ) in the following form 

Ki.H = b\^i,i, + bfii^+bfX,^, (5.17) 

oo 



^hhh = / ^'^^ Mjfef (x)6;f (x)6;j)(x) + {hM.h) perm.} , (5.18) 



2,{ST:fA^nB 



Ml 

= / x'dx M2[bt:^ {x)bf^ {x)b^;i\x) + {hMM) perm.} , (5.19) 



M2 



487r2(nB + 3)(2nB +3)/9W ' 

.AA2{6[7)(x), 

3(87r)3^%B 
487r2(3nB + 3)(2nB + 3)pW 



^hhh = J ^'^^ Af3{b^^\x)b\''J{x)bp^\x) + {hM,h) perm.} , (5.20) 
3(87r)3A34"s+3 

= 48.2(3.. 4)p(^)- ■ ^'-''^ 



The seven-year WMAP data |^] give ~ 2.441 x 10^^. Hence, from (5/?) we can estimate 
the amplitude of the primordial curvature bispectrum F(^{ki, /c2, ^3) is of the order O(10~^^). 
On the other hand, the numerical calculations about the magnetic power spectrum show 
that the amplitude of Pa^ is approximately of the order O(10~^^), i.e. A^^^ ~ 0(10^^). 
Since the some quantities related to the PMFs such as A^, and vr^ have the essential non- 
Gaussian characters, the primordial magnetic bispectrum Fas/ttb (^1 ) ^2, ^3) is proportional 
to the cubic of A^ and vr^ . So, for the magnetic density sourced mode, the amplitude of 
bispectrum FAg{ki, ^2,^3) ~ O(10~^^) is comparable with that of the primordial curvature 
one F(^{ki, k2, k^) ~ O(10~^^). In fact, this observation is just our motivation for this paper. 

5.2 Numerical results 

In this subsection, we will present our numerical results about the reduced bispectra. For 
the case with ns = —2.9 and B\ = 9 nG, the momentum integral (|5.13D is approximately 
equivalent to ( ^.14]) , while ( |5.16| ) diverges in the infrared limit (A; ^ 0). In Figure (^), (|5|) 
and (|^, we plot 6["^(~ b^f^y^op^ and the absolute value of bf^ over / with r]Q = 14.38 Gpc 
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Figure 4: This figure shows 



/(/ + l)&["'(a;)/27r 



10000 



X 10^ for several different comoving radial distances 
X = {jiQ — 0.4?7h.) ~ (t?o ~ 1.677*), where we set the conformal time at present 770 — 14.38 Gpc and at 
the recombination epoch 77* = 284.85 Mpc, respectively. Parameters for PMFs arc 77^ — —2.9 and 
= 9 nG. 



and ?7* = 284.85 Mpc being the conformal time at present and at the recombination epoch, 
respectively. For ub = —2.9, the momentum shapes in b^^^^^ are nearly scale invariant, so 
the profile of -curve (Figure looks similar to those of Ci. However, the main difference 
between 6^"^ and C; is that the former changes the sign, while the latter does not. The 
reason lies in that oc ji{kx)gTi{k), but Ci oc IgTiik)]"^. As what happens to the standard 
L06| , 107], from Figure (^) we can see that the phase of in the high-/ 



model |7g, |To! 

regime oscillates rapidly with respect to x, which will heavily suppress the integrations 
( 5.181) , ( |5.19| ) and (|5.2C1| ) at small scales. However, for the case with the magnetic index 
hb — —3, the integration ( ^.16 ) is dangerous. Because the exponent of the shape k^"-B+^ 
is much less than zero, the full integrand will diverge in the Infrared (IR) limit, if the 
spherical Bessel function ji{kx) and transfer function gTi{k) cannot provide an enough 
positive power to compensate the negative slope of /c3"--b+6^ q^^ numerical calculation 
indeed shows that such divergence does exist. Figure (P) illustrates the IR divergence of 
b\ , which will dominate over all the other terms in the low-/ regime. However, since we 
are interested in the acoustic signatures of bispectrum induced by PMFs, i.e. the moderate 
high-/ regime (/ > 100), we need not worry about the momentum divergence. 

Having the numerical results about the momentum integrations ( 5.13| )- (|5.16D , we can 
finally perform the x-integrations to obtain the reduced bispectra. In Figure (0) and (|8|), 
we plot the integrals 



hih + l)/3(/3 + l)!x^dx b[f{x)b\'^>ix)b\^Jix)/i27T 



X 10^ and 



/i(/i + l)/2(/2 + l)Jx^dx b'^i^\x)bfJ{x)bYJ{x)/{2-nf x 10^ over h by fixing h = h -- 
11, 110, 200, 500, respectively. In the numerical calculations we integrate x from (770 — 2r/^, 
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Figure 5: This figure shows ^[^•'(a;) x lO^'^ for several different comoving radial distances x = 
{r]o — 0.4?7,) ~ (770 — 1.6?7*), where we set the conformal time at present 770 = 14.38 Gpc and at 
the recombination epoch 77* = 284.85 Mpc, respectively. Parameters for PMFs are ns = —2.9 and 
Bx = 9 nG. 



to (?7o ~ 0.1?7*), since the primary signals come from the recombination epoch r/*. From 
Figure (0) and (^), we can see that in the Sachs- Wolfe (SW) regime (/ < 10) our result 
presents a SW plateau, which is consistent with that in [|l^. And in the high-/ regime, 
the integral shown in Figure (^) is greatly damped after the prominent first acoustic peak, 
since the phase of bp\x) oscillates rapidly as a function of x. In Figure the integral 
also has a first acoustic peak, but damps more slowly. And more importantly, from Figure 
(P) we can see clearly that the modes with different parameter sets (Zi, I2) oscillate around 
different asymptotic values because of the effect of the Lorentz force, which is exerted 
by PMFs on the charged baryons. This feature is much different with the one from the 
inflation scenario [0, |10|, |0|, lO^. In order to illustrate this ph enomenon more clearly, 
we perform the above calculations without considering the Lorentz force term in the baryon 
velocity equation ( p.29| ). Figure (^) shows that, contrary to those with the Lorentz force, 
all modes oscillate around zero in the limit of large l^. 

In Figure (|lO|) and (|ll|), we plot the absolute values of integrals ^2(^2 + 1)^3(^3 + 



l)Jx^dx \bij\x)bl^\x)bij\x)\/{27,r\ and [hih+lMh+l] f x'dx |6j;^(x)6j;^(x)6™(x)|/(2vr) 
as a function of I3 with h = h = 11,110,200,500. Since b^'^j^^^ contains the factor b^^\ 
its amplitude experiences a great suppression in the high-/ regime as the same as b^^\ 



From Figure (|10|), we can see that the amplitude of first acoustic peak is approximately of 
the same order as the one in Figure (0), however, in the very high-/ region (/ > 1000) its 
amplitude damps in a power law form. In addition, from Figure (P) and (^), we can see 
that the amplitude of b^"'^ is smaller than that of 5^"^ by the order of O(10~^). Further- 
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Figure 6: This figure shows 1(1 + l)\b^^\x)\/2'!r for several different comoving radial distances 
X = {r]o — 0.4?7,) ~ (?7o — 1.6?7*), where we set the conformal time at present rjo = 14.38 Gpc and at 
the recombination epoch r;* = 284.85 Mpc, respectively. Parameters for PMFs are ub = —2.9 and 
Bx^9 nG. 
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Figure 7: This figure shows the integral \l2{l2 + l)h{l3 + '^) J dx b[f {x)b'i^\x)b'i]\x) / {277^ 
as a function of ^3, with several parameter configurations {li = I2 = 11, 110, 200, 500). 



xlO-^ 



more, for the set of parameters {ub = —2.9, Bx = 9 nG), the coefficients in ( 5.18D -( p^ ) are 
J\fi ~ 3.5 X 10~^^, 7V2 ~ — 6 X 10~^^ and 7V3 ~ — 7 x 10~^^. After considering the hierarchies 
between A/3 and TVi (or A/2), we can safely neglect the ^/.^^jia t^'^-'^ calculations. 
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Figures: This figure shows the integral \h{h + l)l2{l2 + l) J x'^dx b^f {x)b'l^\x)b\^\x)/{2TTy 
as a function of ^3, with several parameter configurations {li = I2 = 11, 110, 200, 500). 
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Figure 9: This figure shows the integral 



xlO^ 



'hih + lMl2 + l)Jx'dx b[f{x)b[^\x)b[]\x)/i27r)'' 
as a function of ^3 in the case without considering the effect of the Lorentz force. From this figure, 
we can see clearly that all modes oscillate around zero in the limit of large I3. 



As a summary, from the above numerical calculations, we can see that the typical 

^2 and ,3/^-^2 

0(10-201 



amplitudes of the reduced bispectra 1%^^^ 



/(27r)2 are of the order 



which is comparable with the NG signals from primary curvature perturba- 



tions 



70, |105| , 106 , 107 1 . In details, b\^} , dominates the total reduced bispectrum with a 
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Figure 10: This figure shows the absolute value of integral hih + 1)'3(^.3 + 

1) / x^dx \b'i^\x)b'l^\x)b'lj\x)\/(2'K)'^ as a function of Z3, with several parameter configurations 
{h ^l2 = 11,110,200,500). 
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Figure 11: This figure show the absolute value of the integral + 1)^2(^2 + 

l)Jx^dx \b[^Jix)b[^^\x)b[^J{x)\/{2n)^ as a function of Z3, with several parameter configurations 
(l^ =/2 = 11,110,200,500). 



(2) 

positive amplitude in the regime / > 100. For b) / , , its amplitude in the high-/ regime. 



is of the same order as the one of b 
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hhh 



but the sign is negative; while it diverges at the 
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large scales. 



6. Conclusion 



Using the full radiation transfer function, we numerically calculated the CMB angular 
bispectrum seeded by the compensated magnetic scalar density modes. For PMFs charac- 
terized by the index ub = —2.9 and mean-field amplitude Bx = 9 nG, CMB bispectrum is 
dominated by two primordial magnetic shapes. For the reduced bispectrum 
which is seeded by the "local-type" shape f^^\k,q,p) (|5.9|), both the profile and amplitude 
look similar to those of the primary CMB anisotropics ]70| , 105, |106| , 107], (see Figure ^, |5|, 



1^ and |8|). However, for different parameter sets {li,l2), such "local-type" estimator ^l^/^/s 
oscillates around different asymptotic values in the high-/3 regime because of the effect of 
the Lorentz force, which is exerted by PMFs on the charged baryons (see Figure I and 
|9|). This feature is different from the one of inflation scenarios where all modes approach 
to zero asymptotically in the high-/ limit. On the other hand, the second magnetic shape 
f^'^\k,q,p) ( |5.10| ) appears only in the primordial magnetic field model. However, the am- 
plitude of the Komatsu-Spergel estimator b^j^^i^i^ ( ^.19 ) sourced by the shape f^'^\k,q,p) 
diverges in the low-/ regime because of the negative slope of shape. In the high-/ regime, this 
amplitude is approximately equal to that of the first estimator b^^}^i,^, but with a reversal 
phase. 

In this paper we only calculated the magnetic angular bispectrum with parameters 
{ub = —2.9, Bx = 9 nG). In fact, the results are strongly dependent on the magnetic 
parameters, especially on the magnetic index ub- Take ub = —2 as an example, the upper 
bound of comoving magnetic mean-field amplitude becomes larger, Bx ^ 25 nG [18|. And 
more importantly, there will appear new Gaussian and non-Gaussian features in the CMB 
anisotropics. At the Gaussian level, the two-point function will have more powers on the 
small scales, i.e. a blue tilt power spectrum; at the non-Gaussian level, the bispectrum will 
no longer diverge in the IR limit. Finally, we would like to comment on the bispectrum from 
the compensated magnetic anisotropic stress mode. Although the amplitude of anisotropic 
stress TT^ff^ is approximately larger than that of density contrast A^^^ by a factor 3, vr^,^^ 
mode only appears in the high order terms in the tight-coupling expansion ( [4.23D -( p:.31| ). 
This results in that the amplitudes of power spectra from tt^^ mode are smaller than those 
from a[.^^ mode, (see Figure |l], ^ and |3|). So we can estimate that the bispectrum from the 
compensated magnetic anisotropic stress should be smaller than the one from the magnetic 
density. 
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A. The covariant approach to cosmological perturbations 

In this Appendix we briefly review the covariant approach to cosmological perturbations 

], H, 1^, H, H, 1^]. Especially, the 



|73, ^, |7^, |76|, |77|, ^ |79|, gO|, |81|, |3|, |8^, 

cosmological covariant formalisms with magnetic fields are carefully studied in |73, 

93, 94, In order to determine the time direction, we define a unit timelike 4- velocity 

vector tangent to the worldline of the observer 



dx°- 



UaU 



-1 , 



(A.l) 



where r is the proper time of the fundamental observer. Then we introduce an orthogonal 
tensor hab = gab + UaU^ with respect to Ua to define the space direction at each spacetime 
point. Using the vector field Ua and projector tensor hab-, we can decompose any space- 
time quantity into its irreducible temporal and spatial parts. Moreover, we can also use 
these fields to define the covariant time and spatial derivatives of any tensor field S^^^,,^'^'" 
according to 



'-'ab--- — " ^ eiJab- - 



(A.2) 



respectively. 

In this paper we use the convention about the effective volume element Eabc and the 
spacetime volume element ??afecd as 



^abc — VabcdU ; 



(A.3) 



where the totally skew pseudotensor is defined as rj^^'^^ = [— det{gab)] ^^"^ ■ Furthermore, 
^abcd jg parallelly transported {jf^'^'^.^ = 0) and satisfies some algebras as 



^abcd^ 

V rjefgh 

„abcs^ 

V Vefgs 



-4! , 



(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 



where the square bracket in the superscript represents for the antisymmetric part of cor- 
responding tensors. 

In General Relativity, the local gravitational interaction is described by Ricci ten- 
sor Rab, while the non-local long-range interaction, such as gravitational waves or tidal 
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forces, is determined by Weyl conformal curvature tensor Cabcd- The decomposition of the 
gravitational field into its local and non-local parts is given by 

Rabcd = Cabcd + -^{gacRbd + 9bdRac — dbcRad — dadRbc) — -^R{gac9bd - 9ad9bc) , (A.9) 

where Weyl tensor shares all the symmetries of Riemann tensor and it is trace-free C^^^^^ = 
0. Furthermore, we can define the irreducible electric and magnetic parts of the Weyl 
tensor 

Eab = Cacbdu'u^ , Hab = -£:''CcdbeU' • (A. 10) 

Then the Weyl tensor can be rewritten by these two tensors 

Cabcd = {9abqp9cdsr - 'r]abqp'r]cdsr)u''u^ - {Vabqpgcdsr + gabqp'ncdsr)^'^^^ , (^.H) 

where gabcd = QacQbd - gadQbc- 

The energy- momentum tensor of a general (imperfect) fluid defined by the observer Ua 
can be decomposed into 

Tab = PUaUb + phab + ^Qi^aUb) + T^ab , (A- 12) 

where p = Tab^u'', p = Tabh''''/3, Qa = -h^Tbcu'^ and Tr^b = h^^%^'^T,d = ^a^^^'^^^d - 
^h'^'^Tcdhab are the energy density, isotropic pressure, energy-flux and anisotropic stress 
tensor of the fluid, respectively. 

In order to characterize the observer's motion we need to decompose the 4-velocity 
gradient into the following irreducible kinematical quantities relative to the Ua-congruence 

VbUa = Tab + ^ab + ^©^afc - ^aUb , (A.13) 

where aab = D^^^Ua), ujab = ^[b^fa]) = Vaw" = Dau'^ and Aa = iia = u'^VbUa are 
respectively the shear and vorticity tensors, the volume expansion scalar, and the 4- 
acceleration vector. In addition, it is useful to deflne the vorticity vector ooa = Sabc'^^ l'^ 
(with {jjab = Sabc^'^) instead of the vorticity tensor. 

A.l Linearized Einstein equations 

Dynamical equations: 

p+ip + p)@ + D''qa = 0, (A.14) 
® + l®^ + \iP + ^P)-D''Aa = 0, (A.15) 

qa + ^@qa + ip+ p)Aa + DaP + D^'Wab = , (A.16) 

2 1 

W(a) + 30w„ + -curl Aa = 0, (A.17) 
2 1 

<^{ab) + -^&<^ab + Eab - ^T^ab - ^(a^6) = , (A. 18) 

1 

2' 

+ ^7i-{a6> + ^^(a?6> + ^^T^ab = , (A.19) 

H^ab) + QHab + curl Eab - ^CUrl TTab = . (A.20) 



E{ab) + QEab - curl Hab + TiiP + P)^ab 
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Constraint equations: 



Daoj" = , 

h 2 

D Gab - curl UJa " g-^a© + = , 
curl (Tab + ^(a^6) " -f^afe = , 
Z)*^,fe + - ^DaP + ^Qqa = , 

D^Hab + ^CUrl ga - (/O + P)uja = . 

A. 2 Two key variables 

It is convenient to define two key variables in the covariant approach 

A« = ^^dJ>^ , Za = aDaQ , 

where 1 = ^,1^, b, c. Taking the spatial gradient of the density evolution equation 
we arrive at 

p(*)A« + (p« +p^'y){Za + aQAa) + aDaD'ql;^ + aQDaP^'^ -p«eA, = , 
For Za, by virtue of the Raychaudhuri equation ( A.15| ) we have 



(A.21) 
(A.22) 
(A.23) 

(A.24) 
(A.25) 



(A.26) 
( WM , 



(A.27) 



• 2G^ 1 , 3 ^ 



ie2 + i(p + 3p) 



Aa - aDaD^Ah = . (A.28) 



B. Equations for matter components 

Under the ideal MHD approximation, the energy-momentum tensors for the five matter 
components are 



T^b = P^'^^aUb + P^'^Kb + 2qf\,) + vrS , (B.l) 



with i = 7, z^, 6, c, and 



l-B'^UaUb + ^B'^hab 

2 



+ vr. 



(B) 

ab 



(B.2) 



Since the total energy-momentum tensor is conserved V^Tab = 0, for each component we 
have 

V'r« = J«=S«n. + M», 5^J«=0, (B.3) 

i 

where E^^^ is the energy transfer and Ma^ the momentum transfer for the i-species. For 
simplicity, in this work we assume the energy transfer vanishes at the linear order (E^^^ ~ 0), 
this gives the energy conservation for each matter component 



+ G(p«+pW) + D"g» = 0. 



(B.4) 
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B.l Photons 



For photons, it is convenient to expand the total intensity brightness I{E, e^) in terms of 
the spherical multipole 



(B.5) 



1=0 



where e'^ is a unit spacelike vector orthogonal to Ua- For CMB, it is usual to define the 
bolometric multipoles by integrating over energy without loss of information 



with 



A, 



47r2'(n)^ 
(2Z + 1)! ■ 



.(7) 



The first three multipoles are respectively 

^ — P ) -'a — Va 1 ^ab — "af, • 

The Boltzmann hierarchies for the total intensity of photons are 



(B.6) 
(B.7) 

(B.8) 



-UeCTT 



(21 + 1)' 
4 1 



(B.9) 



where the right hand side terms stand for the Thompson scattering. The first three mul- 
tipole hierarchy equations are listed as follows. 
Monopole (/ = 0): 



(B.IO) 



Usually one uses the spatial gradient of energy conservation equation, instead of ( [B.10| ), 



A(:') + ^{Za + aeAa) + ^D,D\l 



Dipole (/ = !): 



,(7),(.)_,(7) 



Quadrupole (/ = 2): 



.(7) 



15' 



10 



,(7) 



(B.ll) 



(B.12) 



(B.13) 



Note that the monopole equation ( [B.10| ) is equivalent to the equation of energy conservation 
of the photons, and the dipole one (|B.12 ) to the momentum conservation equation with 
the Thompson scattering. 
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B.2 Massless neutrinos 



For massless neutrinos, the Boltzmann hierarchies are similar with the ones of photons, 
except that the latter has the Thompson collision term, 



4 J, / 

Gai + -^'^Gai + D GbAi + (2/ + + 

4 8 

-GAaj^5ii + — G(Taia2<J«2 = . 

The monopole, dipole and quadrapole equations are, respectively, 
monopole (/ = 0) 



or replaced by 



dipole (/ = 1) 



(B.14) 



(B.15) 



(B.16) 



(B.17) 



quadrupole (1 = 2) 



(B.18) 



B.3 Baryons 

For baryons, the energy conservation equation gives 

and by taking the spatial gradient of ( B.19| ), we obtain 

p('')Ai^) + (p^'') +p('))(Z, + aQAa + aDaD'v'i^) + aQDaP^''^ - p^'^QA^^^ = . (B.20) 



The momentum conservation reads 



(B.21) 



where the momentum transfer for baryons is due to the Thompson scattering and Lorentz 
force from PMFs. 
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B.4 Cold Dark Matter 



Since CDM only gravitates, the density contrast equation is given by 



(B.22) 



or 



(B.23) 



and the velocity equation is 



(B.24) 



They are simple. 
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